Stellar Relaxation Time 


[Chandrasekhar 1960, Principles of Stellar Dynamics, Chap II] 
[Ostriker & Davidson 1968, Ap.J., 151, 679] 


Do stars ever collide? Are interactions between stars (as opposed 
to the general system potential) important? We can answer this 
question by calculating the time it takes for a star’s orbit to 
be “significantly” perturbed by individual encounters with other 
stars. To calculate this relaxation time, let’s first define the word 
“significant”. One way of doing this is through total energy: 
when does the kinetic energy exchanged during stellar encounters 
equal the star’s original kinetic energy, t.e., 


Ty => ) (AE) =E (30.01) 





But for simplicity, we’ll define “significant” as the time it takes 
a star to lose all memory of its original trajectory, i.e., 


Tp => J osin?y=1 (30.02) 


We then assume that a) all deflections are two-body encounters, 
b) each encounter is statistically independent, and c) close en- 
counters are insignificant compared to long-range encounters, so 
that during each encounter, |AE| « E. Under these assump- 
tions, all the deflections are small (siny < 1), and we can use 
the Born approximation, where (Vinit © Ufnal © V). 
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For a single encounter, the deflection angle, y is related to the 
initial impact parameter, b, by 
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From the geometry of the encounter 


F, = Fsin0 = F (2) = (<a @ (30.04) 


Also, from the Born approximation 
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So 
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or, since r = (s? + b?)1/2 
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Letting xz = s/b 
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Since for small deflections, tan Y © y = v/v 


y= (30.09) 


Now, let’s sum this over all possible collisions. The number of 
collisions that take place in time dt depends on the impact pa- 
rameter, the distance a star travels in dt, and the density of stars 
in the stellar system, N, i.e., 


Neoll = (27b db) : (vdt) -N (30.10) 


So, to deflect the star by 90°, 
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As for the limits on the log quantity, we can use the obvious fact 
that no deflection angle can be greater than m. Thus 
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Similarly, it is clear that the maximum impact angle must be less 
than the mean distance between stars, so 
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Obviously, the above derivation involves a number of approxima- 
tions. A more rigorous derivation by Chandrasekhar gives 
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where H(x) and G(x) are factors of the order unity that depend 
on the stellar distribution function. Finally, Ostriker & Davidson 
(1968) give an improved, recursive expression for the relaxation 
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In the solar neighborhood, this timescale is much larger than a 
Hubble time. Thus, the motions of stars are “collisionless”, and 
controlled only by the overall Galactic potential. 


The Collisionless Boltzmann Equation 
[Binney & Tremaine 1987, Galactic Dynamics, 1987] 


The basis of understanding galactic dynamics lies with the col- 
lisionless Boltzmann equation. Imagine a closed volume, V, 
bounded by a surface $, and containing a mass, M(t). The net 
amount of mass flowing through a surface is equal to the change 
of mass in the volume, i.e., 


[ovas=- —_—_ -| oe dV (30.18) 


But the divergence theorem in mathematics says 





[eas= | V - Qdv (30.19) 
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Now expand this concept to stars flowing in a 6-dimensional 
phase space (2,Y,Z,Uz,Vy,Uz). Assume that the flow of stars 
is smooth, and is controlled by a potential per unit mass ®. Now 


let 
f = the stellar density in the 6-D space 


w = the 6-D position coordinate 


x = the 3-D space coordinate 
v = the 3-D velocity coordinate 


With these definitions, the continuity equation for a fluid of stars 
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Now let’s expand this out to 


ay Coron aa i | =0 (30.23) 


V- (fw) =0 (30.22) 


and explicitly write out the space and velocity parts of the second 
part of the equation. 
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Since z and y are independent dimensions, 
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Also, for conservative forces with no collisions 
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Also, since iie ee a ae feels depends only on its lo- 
cation, and not on its velocity 
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Finally, to make the equation more transparent, we can again 
explicitly write out the eas and a parts of the sum 


ofa rage 530 st =0 (30.30) 














or, if we substitute t = v and ù = Toe 
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This is sometimes called the collisionless Boltzmann equation, 
the Vlasov equation, or the equation of continuity. Note that 
this equation, when expressed in terms of a Lagrangian derivative 
(i.e., from the point of view of an observer moving through space 
with the fluid) is 


Ov Ov df 
w= (Z) a+ (Z) a = Pree (30.33) 





Jeans Equations 


As it stands, the collisionless Boltzmann equation is rather use- 
less, as it is not only a function of 7 variables (many of which are 
not easily observable), but their derivatives. Even if f(wW) were 
measurable, uncertainties due to Poisson noise would play havoc 
with the derivatives. Fortunately, more tractable equations can 
be found by finding the moments (averages) of the equation. Such 
moments produce the Jeans Equations. 


First, note that the spatial average of any quantity associated 
with the 6-D stellar fluid is simply found by integrating over 


velocity, i.e., 
Q) - [ores] fren (30.34) 


Moreover, since the density of stars in the fluid, v, is simply 
v(z) = J fdg (30.35) 
the spatial mean of any quantity is 
1 gu 
=— | Q- fdv (30.36) 
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So, for the first Jeans equation, start with the Boltzmann equa- 


tion 
Ob əf 
ory DT o 5 | = 0 (30.31) 


and integrate over the 3 velocity dimensions 


Fy ofeg dy df Ze v=0 (30.37) 

















Since time, position, and velocity are independent quantities, we 
can extract those dependencies from the integrals 
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The first integral in the equation is simply the density, and the 
second is just the density times mean velocity, so (30.38) becomes 
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But if you expand that last integral 
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(because obviously, the phase-space density of objects at infinite 
velocity is zero). This leaves us with the first Jean’s equation: 
the 3-D equation of continuity 
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For the second Jeans equation, we again start with the Boltz- 
mann equation 
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This time, however, we multiply through by v; before integrating 
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(30.42) 
Once again, the time and space derivatives can be taken outside 
the integral, leaving mean quantities 
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And, once again, we can expand and evaluate the individual in- 
tegrals of the last term 


fo By = JII v SL duidvjdey (30.44) 


If 7 Æ i, the result is the same as before 


[fos in dv;idvjdvk = Jj vjdojdos f df =0 (30.45) 
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But if 7 =2, 


JIJ viS duidujdus = Jj dvjdey | vidf (30.46) 


We can evaluate this last term by integrating by parts 
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The first term is zero, and the second is just the space density, 
v, so the second Jean’s equation becomes 
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The final Jeans equation comes from noting that the anisotropic 
pressure term (i.e., the stress tensor), oj; is 


Gij = ((vi — (vi) (vz — (vj))) = (very) — (vi) tvz) (30.49) 
So, if we substitute this into the 2nd Jeans second, 
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If we expand this out 
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multiply the 3-D continuity equation by (v;) 
ðv wa, ô (vv) 
D ez 2 p = 0 (30.52) 


and subtract the two equations, we get 
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Here, vo;; is the anisotropic pressure (stress) tensor. But since 
Gij is symmetric, its matrix can be diagonalized. This produces 


the principle axes of the velocity ellipsoid. 


This is the equivalent of Euler’s equation for fluid flows 


Ov 


Pap +p(v-V)v = pV®— Vp (30.54) 


where p is the density and p the pressure. 
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Jeans Equations in Cylindrical Coordinates 


The Jeans and Boltzmann equation are almost never used in 
their Cartesian coordinate form. For most applications in spiral 
galaxies, cylindrical coordinates are used; in elliptical galaxies, 
the equations are given in spherical coordinates. The derivation 
of these equations is straightforward, but tedious. For the cylin- 
drical coordinate equations, the Boltzmann equation becomes 
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and if we assume azimuthal symmetry, the Jeans equations are 
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This last equation is often used in its simplified form. First, 
assume the galaxy is in a steady state, so that the time derivative 
is zero. Next, consider (vgv,). From symmetry, this term should 
be zero in the plane of the galaxy. Above and below the plane, 
this term might be zero, but worst case is that the principal 
axes of the stellar velocity ellipsoid are rotated to align with that 
of the spheroidal component. To estimate the effect this would 
have on (vgu-z), we must first transform the spherical coordinate 
system (Vr, Vg, Vg) to the cylindrical system (vp, vg, vz). This is 
tedious, but straightforward. If the principle axis of the system 
is aligned with the disk (i.e., if (v-vg) = 0), the the result is 


(uRvz) = { (up) — (v9) } Roe 


9 9 2 1 
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Since most of the stars are near the plane of the galaxy, uv, ~ UR 
and vg & vz, SO 


(uve) © {(u2) — (v2)} - = (30.60) 


Now, when we evaluate the third term of the Jeans equation 
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Since most of the stars in a disk galaxy are near the galactic 
plane, i.e., z < R, this term is small. So except for the region 
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near the galactic center 





~ 0 (30.62) 
Similarly, the derivative 
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(And remember — this is the worst case scenario. The closer to the 
plane you are, the more likely that (v,v,) = 0 by symmetry.) So, 
to first order, the Jeans equation relates the z velocity dispersion 
to the galactic potential by 
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To see the importance of this, we can take the derivative of this 
equation, 
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and then note that near the galactic plane, the gradient of the 
potential is almost entirely in the z direction. Hence 
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via Poisson’s equation. Thus the mass in the disk can be mea- 
sured via the z-motions of stars near the galactic plane. 
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Jeans Equations in Spherical Coordinates 


The Boltzmann equation in spherical coordinates is 
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and, for steady-state systems with (v,) = (vg) = 0, the Jeans 


equation is 


PEE A Y faw) — (pte) (30.68) 


Another way to look at this equation is to think of it in terms of 
hydrostatic equilibrium. Recall that the definition of pressure is 


P = zp?) = plv?) (30.69) 
where v; represents one component of the motion. So let P be 


the radial pressure imparted by the stellar motions, and Q be 
the tangential pressure term, t.e., 


P = vlo?) Q =v ((v5) + (va?) (30.70) 





In this case, the Jeans equation looks is 
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If the stellar orbits are randomly distributed (i.e., isotropic, so 
that ((uz) = (vg) = (v3), then 2P = Q, and the equation re- 
duces to the simple hydrostatic equilibrium. 





To describe the degree of orbital anisotropy in a spherical system, 
one often uses the parameter $, 
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For isotropic orbits, 6 = 0; for purely radial orbits, 6 = 1, and 
for purely circular orbits, 6 = co. Note that 6 need not be 
a constant throughout the system; realistic models of spherical 
systems often have {(r) as their most important free parameter. 
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